Evaluation of pairwise entanglement in translationally invariant systems with the 

random phase approximation 
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We discuss a general mean field plus random phase approximation (RPA) for describing compos- 
ite systems at zero and finite temperature. We analyze in particular its implementation in finite 
systems invariant under translations, where for uniform mean fields it requires just the solution of 
simple local-type RPA equations. As test and application, we use the method for evaluating the 
entanglement between two spins in cyclic spin 1/2 chains with both long and short range anisotropic 
Xy-type couplings in a uniform transverse magnetic field. The approach is shown to provide an 
accurate analytic description of the concurrence for strong fields, for any coupling range, pair sep- 
aration or chain size, where it predicts an entanglement range which can be at most twice that of 
interaction. It also correctly predicts the existence of a separability field together with full entan- 
glement range in its vicinity. The general accuracy of the approach improves as the range of the 
interaction increases. 

PACS numbers: 03.65.Ud, 03.67.Mn, 75.10.Jm 



I. INTRODUCTION 

The random phase approximation (RPA) [H-ll] is a 
well-known technique in many-body physics. R can be 
considered as the next step after the mean field approxi- 
mation (MFA), being able to describe in a rather simple 
way some of the efi'ects induced by the residual interac- 
tion, such as collective excitations In this work we 
want to examine its application to the problem of eval- 
uating pairwise type entanglement in general composite 
systems invariant under translations, such as cyclic spin 
chains with long or short range couplings in a uniform 
magnetic field, at both zero and finite temperature. The 
fundamental importance of quantum entanglement in dif- 
ferent areas of physics is well recognized, constitutiiig an 
essential resource for quantum information science [j, Q 
and providing a deeper understanding of quantum cor- 
relations in many-body and condensed matter physics 
0-i8i]. Nonetheless, the evaluation or even the estimation 
of entanglement in interacting many-body systems is in 
general not an easy task, particularly for long range cou- 
plings and finite temperatures, lying beyond the scope of 
basic methods like the MFA which rely on separable trial 
states. 

Here we will show that the MF+RPA can provide a 
simple general method for estimating pairwise entangle- 
ment, with a complexity which does not exceed that of 
solving a local MF+RPA problem in the case of transla- 
tionally invariant systems with uniform mean fields. Rs 
accuracy actually increases for long range interactions 
or high connectivity, i.e. for situations where numeri- 
cal techniques for evaluating ground states of spin chains 
(like Quantum Monte-Carlo (9,], DMRG ^] and meth- 
ods based on matrix product states [ll|) become nor- 
mally more complex to apply or less accurate. In any 
case it allows for a rapid estimation of the main features 
and their behavior with the control parameters, leaving 
the application of more accurate approaches for a second 



step. We have previously shown that for fully and sym- 
metrically connected spin systems (Lipkin-type models 
[H [H), a MF-hRPA treatment is indeed able to de- 
scribe the pairwise entanglement at both zero and finite 
temperatures [l^ | , becoming exact in the thermodynamic 
limit. Its ability to reproduce pairwise entanglement in 
more general systems was, however, not examined. 

We will first briefiy revisit the general MF-I-RPA for- 
malism derived from the path-integral representation 
of the partition function, discussing its implementa- 
tion in composite systems and in particular in those 
which are translationally invariant. We next apply the 
method to finite cyclic spin 1/2 chains with general 
range anisotropic XYZ type couplings. Comparison 
with numerical exact results is made for finite chains 
with anisotropic XY interactions of distinct ranges. The 
method is able to capture most essential features of the 
entanglement between two arbitrary spins away from MF 
critical regions, becoming accurate for strong magnetic 
fields, where it provides an analytic description of the 
concurrence. At weak fields the agreement with exact 
results is less accurate but improves as the interaction 
range increases, being as well able to predict the appear- 
ance of a factorizing field and an infinite entan- 
glement range in its vicinity. 



II. FORMALISM 
A. General Mean Field+RPA treatment 

We consider a general system of n distinguishable con- 
stituents with Hilbert space dimensions d,, interacting 
through a general quadratic Hamiltonian 



(1) 



where we have adopted tensor sum convention for re- 
peated labels and stand for general independent lin- 



2 



ear combinations of ZocaZ operators, i.e., = Yl7=i 



with 0„j = Ii ( 



H+l 



ii i ^ j). We will assume V^" = V^^, as commutators 
[Ofj,, Ou] are again linear combinations of local operators 
and can be included in the linear term in ([ij . A hamilto- 
nian linear in represents obviously a non interacting 
system, being diagonal in a basis of separable states and 
requiring just di x di local diagonalizations, whereas H 
demands in principle a x di diagonalization, its 

eigenstates being entangled in general. 

The partition function Z = Tt exp[^/3H] admits, how- 
ever, an exact representation in terms of linear hamilto- 
nians by means of the auxiliary field path integral 



Z = 2?[0]Trf exp 



HmT)]dT , (2) 



(3) 



where (j)'^ are the auxiliary fields and T denotes (imag- 
inary) time ordering. The operator under the trace 
in Eq. ^ is just the imaginary time evolution opera- 
tor associated with the linear Hamiltonian H[(I){t)], be- 
ing then a product of local evolution operators. Eq. 
^ can in principle be evaluated through a Fourier ex- 
pansion 0'^(t) = E™=-ooCe''"""/'', with V[cb] = 
n™[det(^)-V2n^d0M]. 

The MF-hRPA treatment, to be abbreviated as CMF 
, is obtained by evaluating Eq. 1^ in the gaussian ap- 
proximation Q around the static mean field 4>l^ — 5mo4'^ 
which maximizes 



Zmf(0) = Trexp[-/3iJ(0)] = 6^2 



where Zi((j)) — tr exp[— /3(5'^ — 0'')opi] is a local partition 
function. It then satisfies the self-consistent equations 



9 In Zi{(f)) 



(5) 



such that {H((j)))^ = bf'iO,,)^ - ^V^''{Of,)^{0,)^ at a 
solution. The final result can be written as 

ZcMF = ^mf(0)C'((/)) , (6) 

OO 

= n Det[5^; + (7) 



i TT ! 



u!a sinh 

a>0 sinh ' 2 



(8) 



where i?™^, = X]j=i ''"^wU) MF response matrices and 
2nim , , , , , , 



(9) 



fiP) 



are local responses {r^^{j) = —d{Of_ij)^/d(t)'' is a local 
susceptibility matrix), with {b'^ — 4>^)o^j_j\K,j) — eK-\Kj), 

= S^.^', and = e'^"''^ / Zjicj)). C{(f)) contains 

the static (<50o) and quantum ((5?!'^_^g) gaussian fluctua- 
tions around the mean field and requires just local diag- 
onalizations if evaluated through Eq. ([7]). In the closed 
form ([5]) , Aq = £„j — e^'. , with a labelling all pairs of dis- 
tinct local eigenstates (a > indicating Kj > k^), while 
oja are the RPA energies, obtained from 



(11) 



where i?(w) = Sj^Oj'^)- They are the poles of the 
RPA response matrix [/ -I- R{uj)V]^^R{u!) and come in 
pairs of opposite sign. They can also be obtained as the 
eigenvalues of the RPA matrix 



(12) 



where /« = p^, . - 0^^ = ('^i lo/xj l'*j>7 of dimen- 
sion dj{dj — 1). Let us mention that under a linear 



transformation On 



V^O^, we have b^' 



yt^u _^ yt^u ^ u^U^VP^, Eqs. dZ])-® being of course 
independent of the representation. 

Eq. ^ can be applied away from MF critical points 
(where the static determinant in (O-® and the lowest 
RPA energy will vanish, and where the approach can be 
improved for T > by integ rating exactly over the rele- 
vant static variables [3, iUll^l), becoming accurate for 
small VR. In the presence of vanishing RPA energies 
arising due to a mean field which breaks a continuous 
symmetry of H ^3|, the product in ([8]) remains finite but 
/i, ly in should be restricted to the intrinsic static 
fields, with static orientation variables integrated out ex- 
actly and contributing with a prefactor to ([5]) If 
7^ V a, we may rewrite ([8]) as 



C(0) = Det[ 



V^R 



■yj sinh 



(13) 



where = [R^ - R{0)][1 + VR{0)]-^ vanishes for T ^ 
and the last factor is just the ratio of partition functions 
of independent bosons of energies Wq and Aq. For T — > 
the energy Eqmf = _d\nZ^ME. approaches the usual 

MF+RPA expression @] (7?(0))0 + \ EJ^o. " ^c)- 

For a Hamiltonian representation in terms of purely 
local operators O^i , we should replace and cj)^ by O^i 
and (j)P^ in previous expressions. 



(14) 



and = <5.,r;r.(*), such that F'^^i?™ ^ l^^^«C.(j). 

Eqs. ©-([HI) will then involve in general determinants 
of matrices connecting all components. We can assume 
ytiiiyi — 0, as self-energy terms are local operators and 
can in principle be also included in the linear term. 

Although the representation (ITi|) is not necessarily 
the most convenient one for evaluating C(0), it allows 
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to evaluate two site averages directly as {0^iO,yj) = 
2/3-^9 In Z/ay'^*''^ leading in CMF to 

(Om.O.,> - {0,.),{0..j), + -^^^^ ■ (15) 

The reduced density matrix for the i — j subsystem can 
then be recovered by considering a complete set of lo- 
cal operators. For degenerate symmetry breaking mean 
fields Eq. (|15p should be averaged in principle over the 
different solutions. 



B. Translationally invariant systems 



or equivalently, the eigenvalues of the effective local RPA 
matrix aaa'{k) = \a5aa' + faOf_,.-aV^''{k)o^a,, of di- 
mension d{d — 1). C{(j)) reduces then to the prod- 
uct of n single site correction factors with couplings 
v^^^{k). These results also hold for D- dimensional 
cyclic systems (for instance spins in a torus) provided 
yt^ii^j _ -y^i^^j _ replacing matrices v{k) by v{k) ~ 

Eq. dTSl) will now depend just on the separation i — j, 
becoming 

/n n \ - \2 I ^ \^ i27Tkj/n9}^l£W ,r,,^ 

{0,.0.,+,)-{o,)^ + — }_^e -g-;^- (21) 



Let us now consider the case of identical components, 
i.e., identical Hilbert spaces {di = d) and operators 
(o^i — o^) at each site, with 5'*' = and 



V 



(16) 



where v^'^{n — j) = v^'^{—j) for a finite cyclic chain. In 
this situation we may conveniently rewrite Eq. (jl4p as 



(17) 



k=0 



where v{k) is the (discrete) Fourier transform of v{j), 

n-l 



j=0 



(18) 



and similarly, O^^ = n X]J=i ^^^'^'^"'^"^w (such that 

Om = Efce-'2^'=^'/"0^fc)- Thus, y^'^'^'^' = n5k,-k'i^''{k) 
in Fourier representation. 

We will also assume a uniform mean field 0^' ~ 0^, 
such that {Ofj_i)^ — {o^)^ and hence (Eq. [S|), 



(19) 



which is an effective local MF equation depending just 
on the total coupling ■C'"^(0) = Notice that in 

Fourier space Eqs. © become 0''*^ = ni^''' {k){6^^^k)4>, 
the uniform solution corresponding to 0^'^ = nS^'~'(j)'^ and 
leading to (iJ(0))0 = n[b^{o^}^ - iw^'^(O)(o^)0(o,)0]. 
In this case r'^^{i) = r™^ is site independent, imply- 

and therefore, {yK^YX' = ^I'^^^^kY^^, diagonal m k. 
Hence, Eq. ([8]) becomes 



c(0) =n 



n 1 -i-r LOnik) sinh 
Dct[5^; + S^''(fc)r°J-^n 



Q>o smn 2 



with Wq, (A;) the roots of the local RPA equation 
Det[Sii + i^P{k)rpAoj)]=0, 



(20) 



III. APPLICATION 

A. Finite spin 1/2 chain with general range 
XY Z-type couplings 



We now consider a finite spin 1/2 cyclic chain in a 
uniform magnetic field. The local operators are the 
spin components s^, n = x,y,z, and we will assume 
'^^^{i) — (5^'^w^(j) (j-independent principal axes), with 
the magnetic field parallel to one of these axes (z-axis). 
This wide class of systems comprises well-known mod- 
els such as the Ising and 1-D XY models with nearest 
neighbor couplings [y, [2l| , as well as the Lipkin model 
0, [l2| - [l^ . where every pair is identically coupled. The 
Hamiltonian reads 



(22) 



where we will assume w''(j) = v^{n — j) — v^{—j). 
Eq. (1^^ always commutes with the "5*2 parity" Pz — 
ll^e'-(S.'+'/^), entailing (S^,) = 0, {S^^S,,) = for 
fi = x,y and j ^ i at any T > 0. 

Eqs. p9)) for a uniform MF become here 



=w^^^tanhi/3A, 



where -D^ ^^(j) 



A ""^"^ 2 

g-i27rfej7' 



(23) 



' - 5^)2 and 

&^ — (0, 0, 6). We will focus on the anisotropic attractive 
> 0, where the lowest solution cor- 
(normal solution). 



case Vq > \v, 



bo 4 


-\b\ 


ba- 


-\b\ 



responds to 0^ = and i) (/)^ 
valid for \b\ > be — Vq — Vq or 

if |6| < be, where (j)^ — — -Dq tanh i/?A, or otherwise ii) 
(j)^ = 7^ {degenerate parity breaking solution), 

where X — v'q tanh ^/3A, (j)^ — —v^b/bc- The ensuing 
CMF treatment involves here just 2x2 diagonalizations 
with a single RPA energy for each value of k: 



,iio^{«M>|(2coshi 



k Z 



(24) 



4 



where, defining 7^, 



A 



and / = tanii hf3X, 



wfc = A^(l - fijl/X)[l - fi^lil + ilvD/X] , (25) 
r:9. - ^ . (26) 



The spin correlation 
ated as 



can be evalu- 



1 91nZ, 



CMF 



n/3 dv^^{i) 



(27) 



where (s,)^ = 17./ and a^^ = ^ e'^-'^.V-^lgp. 
The reduced two-site density matrix in the standard basis 
can then be recovered as 






ttxj — OL. 



yj 



Otxj + Oiyj 












p7 



where = i + a^j ± (sz) and {sz) — —(i~^d\nZ/db. 

We are here interested on the pairwise concurrence Cj 
[2^ , a measure of the entanglement between spins i and 
i+j, given in this system by Cj = Ma,x[Cf , CJ , 0], where 



Ct = % 



ctfr.i - a. 



a 



1/4] 



(28) 



CJ - 2[|a,, + a,, I - ^(l/4 + a.,)2-(s,)2] , (29) 

represent a parallel or antiparallel concurrence respec- 
tively 16]. Just one of C^ can be positive in a given 
state. Note that Cf < at the MF level (a^j = 



for an anisotropic interaction of range L (w±(j) = for 
j > L and v_{L) ^ 0) the T = entanglement range for 
\b\ ^ be can be at most twice the interaction range. Com- 
parison with exact perturbation theory indicates that for 
high fields, Eq. (I3ip is actually exact for any n but up 
to the first non-zero order. For instance, in the near- 
est neighbor XY case v'^{j) — v^^{6ji + Sj^n~i)/2, with 



0, Eq. ^ leads to C+ = if j > 2 and 



Ct 



2b 



C 



\v-\i\v+\-\v-\) 
4&2 



(32) 



with v± — (w^±t;^)/2, which coincide, up to 0{v/b) and 
0{v/b)^ respectively, with the exact result for the concur- 
rence obtained with the Jordan- Wigner transformation. 
Hence, in this limit there will be 0(v/b)'^ concurrence 
between second neighbors if Jf+J > \v-\. 

For r > 0, the main thermal corrections to ([5T|) will 
arise from the decrease of the MF contribution (sp)^ to 
a^j, leading to 



C+(T) « C+(0) 



2e" 



-/3A 



(33) 



for sufficiently low temperatures such that C^iT) > 0, 
where A = 6 + and C+(0) is the T = value (I3I|). We 
have neglected in (j33)) thermal corrections to ajj^ , which 
will lead to higher order terms in v/X. From p3p we may 
estimate the limit temperature for pairwise concurrence 
at high fields. 



A/ln[2/C+(0)] 



(34) 



which will increase almost linearly with increasing b 



B. Concurrence for strong fields 



C. Separability field 



Let us first examine the concurrence for strong fields 
b 3> 6c, where the mean field solution is always normal 
(and uniform) and the ground state is the fully aligned 
state plus small corrections. Eq. becomes 



UJk - A^i - /— - 2/ 

1 lo,^ . 



3 (^-fc 



+ 0(/f)4] 
(30) 



where A = b+fv^, = ^{i'k^^k)^ and we have assumed 
v'^{j) — 0{v). In this regime the exact GS concurrence 
can only be parallel. Up to 0(u/A)^ Eqs. (I27l)-(|30l) then 
lead at T = to 



C^ 



A 



A2 



En — 1 2 / ■ 



2A2 



(31) 

where v±{j) — (w^(j) ±z)^(j))/2. Hence, pairs connected 
by v-{j) will exhibit in this limit a parallel concurrence 
of first order in v/X, whereas those unconnected may still 
exhibit a parallel concurrence of second order in v/X if 
linked by the convolution of v^ with v^ . This entails that 



Let us now assume a common range such that w^(j) = 
r{j)v^, with — 1 (^0 ~ ^'')- Anisotropic chains 

with v^ < v^ < ti^ and r{j) > will exhibit a factor- 
izing field [1B4I3 bs = \J {v^ — v^){vy — v^) < be where 
the degenerate parity breaking MF states become exact 
ground states and Cj vanishes for large n fv^, changing 
from antiparallel {\b\ < bg) to parallel {\b\ > bg). It is 
verified that at T = and b — bg, a^^- = for 7^ in 
p7|. entaihng Cf = V j > also in CMF. Expansion 
of LUk around 6s actually leads at T = to 

'•y b,b~b, ^,b-b 



UJk = u^[l - rfc — 



Tk- 



o{- 



s \2 



where = e ^^^^''r{j), implying, up to 0(^''= 



±7, 



b,b — b„ 



-E- 

n ^ 1 



be v^ 



Tk 



rkvV /v'-^ 



(35) 
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FIG. 1. (Color online) Comparison between exact (solid lines) 
and CMF (dashed lines) results for the ground state concur- 
rence Cj of spin pairs with separation j as a function of the 
transverse magnetic field b, for a long range XY coupling 
v'^{i — j) cx v'^ /\i — j\" , with x = Vy/vx = \ and n = 18 spins. 
The concurrences vanish at the factorizing field 6s = y/X^c 
where 6c = denotes the MF critical field. The inset depicts 
the reentry at high fields of the concurrence of distant pairs 
for a = 2. 



where r™(j) = J2i^U ~ i)r^"'^^{i) {m > 2) denotes the 
TO*'' convolution of r{j). For any finite coupling range 
satisfying r{j) > for 1 < j < L and otherwise, Eq. 
(pel) yields 7^ > for j = 1, . . . , n. Therefore, CMF will 
predict in this case full entanglement range in the imme- 
diate vicinity of bs , with Cj changing from antiparallel 
to parallel as b crosse s 6 c 1 which is in agreement with the 
general exact result The slope of Cf{b) at & = 6^ 

is, however, not necessarily exact in CMF. 



D. Comparison with exact results in finite chains 

Illustrative results for anisotropic XY couplings 
(v^{j) = 0) with different ranges are shown in Figs. [TJ2] 
as a function of the transverse field. We first consider in 
Fig. [Da long range coupling of the form v^{j) cx u''/|j|" 
for 1 < IjI < n/2, with = > 0, where exact ground 
state results for n = 18 spins have been obtained by 
direct diagonalization. We have selected an anisotropy 
X = v^/v^ — 1/2, in which case the factorizing field is 
bs = y/xbc ~ 0.716c, with bc = v^. As predicted by CMF, 
at 6 = 6s the exact concurrence is seen to vanish for all 
a, reaching always full range in its vicinity (for finite n 
the exact result actually approaches at T = exponen- 
tially small a and j-independent finite lateral limits [l7| 

C± = (1 - x) ^^l"/2 at 6 = 6,, with C± « 0.002 for 
n = 18 and x = 1/2, not predicted by CMF). 

The a = case corresponds to the Lipkin model [l3| . 
where v'^{j) = /{n ~ 1) and = v'^{n5ko — — 1)- 
In this case Cf ^C^\f j, with <2/n^ due to the 



monogamy property [2J| . CMF is here quite accurate for 
all fields values away from 6^, providing the exact result 
for the rescaled concurrence nC for large n [3]. 

As a increases, CMF remains accurate for high fields 
6 > 1.56c, where the concurrence is correctly described 
by Eq. (pij) . i.e., Cj oc (w_/6)/|j|". For sufficiently large 
a Eq. pip actually predicts a weak reentry of the con- 
currence Cj at strong fields for large separations j, since 
the last second order term in pTjl will be negative and 
greater than the ffist order term for not too strong fields 
if J is sufficiently large. This reentry is confirmed in the 
exact results for large separations, as seen here for a = 2 
(inset of bottom right panel) . CMF looses precision for 
low fields |6| < 6c, although for a < 1 it is still quite 
reliable for |6| < 6^, where its accuracy increases as j in- 
creases. Notice also that for a < 1 we obtain for n — 18 
full range concurrence at all fields, whereas for a = 2 
the concurrence becomes very short ranged at low fields 
(j ^ 3), being non-zero for large j just in the vicinity of 
6s or at very strong fields, i.e., where the nearest neigh- 
bor concurrence becomes small, in agreement with the 
monogamy property. This behavior is qualitatively re- 
produced in CMF. Let us finally mention that for a — 2, 
results for the first few Cj will remain stable as n in- 
creases (as 1/j" is in this case convergent), those of 
CMF remaining close to those depicted for n — 18. 

Fig. [5] depicts results for finite range couplings of con- 
stant strength, i.e., v'^{j) = ^v^/L for \j\ < L and 
otherwise (such that = v^), at the same anisotropy. 
For nearest neighbor coupling, which corresponds to the 
a — 00 limit of the previous case, exact results for any 
finite n and T can be obtained with the Jordan- Wigner 
transformation [2l| plus parity projection CMF is 

again confirmed to be accurate for high fields for both 
j = 1 and j = 2 (Eq. while for |6| < 6c it pro- 

vides only a qualitative agreement (with correct predic- 
tions like the full entanglement range in the vicinity of 
6s), even though it is still reliable for standard observ- 
ables like the spin correlation axi away from 6c (inset 
in the upper left panel). The thermal behavior of Cj is 
also correctly described by CMF away from 6c, as seen in 
the upper right panel, where exact results confirm the in- 
crease in the limit temperatures Ti and T2 for high fields 
as predicted by Eq. 

Nevertheless, the accuracy of CMF at low fields im- 
proves as soon as the range L is increased, i.e., as 
v^{j)/v^ decreases. For instance, results at 6/6c = 0.1 
significantly improve already for L > 2, as seen in the 
bottom right panel, while for L = 3 CMF is seen to 
provide the correct general picture except in the vicinity 
of 6c (bottom left panel). In particular, the concurrence 
range for high fields is seen to be again twice the coupling 
range, in agreement with Eq. (|3ip (actually, for j — 6 
both the first and second order terms in Eq. PT|) vanish 
for X = 1/2 and L = 3, and an expansion up to 0{v^ /b)^ 
is required, Ce{b) being still positive in both CMF and 
the exact results). The splitting of the concurrences Cj 
for j = 1, 2, 3 is as well a second order effect. 



6 




1 b/bc 2 0.2 T/bc 0.4 




1 2 2 4 6 8 



FIG. 2. (Color online) Results for finite range XY couplings 
at the same anisotropy x = 1/2- Top panels depict the con- 
currence for nearest neighbor coupling (L = 1) and n — 100 
spins, as a function of the transverse field at T = (left) 
and as a function of temperature at fixed fields (right). The 
inset depicts the spin correlation {SxiSx2) at T = 0. Bot- 
tom: Results for interaction range L and constant strength 
{v'^{i~j) oc f'' if |i — j| < Z/ and otherwise) for n = 18 spins 
at T = 0. The left panel corresponds to L = 3 (third neighbor 
coupling) while the right panel depicts the concurrences Cj 
as a function of the range L at fixed low field. 



IV. CONCLUSIONS 

We have examined a general MF-I-RPA treatment for 
describing composite systems with quadratic interactions 
at both zero and finite temperature, showing that it be- 
comes particularly simple for finite translationally invari- 
ant systems with uniform mean fields. The approach is 
capable of reproducing the main features of the pairwise 
entanglement, for all pair separations, in cyclic spin 1/2 
chains with anisotropic XY couplings of different ranges, 
away from MF critical regions. It also provides the cor- 
rect asymptotic behavior of the concurrence for strong 
fields, where it predicts interesting features like the pos- 
sibility of a reentry of the pairwise concurrence for large 
separations, as well as an entanglement range which can 
be at most twice that of the interaction for finite range 
couplings, which were confirmed in the exact results. It 
also predicts the factorizing field and the full entangle- 
ment range in its immediate vicinity. 

The method is specially suited for treating systems 
with high connectivity or long range interactions, where 
its accuracy improves. Let us remark that the individ- 
ual components are in principle arbitrary in the present 
formalism. They could be also chosen as small arrays of 
coupled spins or subsystems treated exactly, leaving the 
RPA for the remaining interactions, a possibility which 
is currently under investigation and which could improve 
results for finite range couplings or dimer type chains. 
The extension to higher dimensions is as well straightfor- 
ward. 
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